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The peculiar motion of the Earth causes a dipole anisotropy modulation in the distant galaxy 
distribution due to the aberration effect. However, the amplitude and angular direction of the effect 
is not necessarily the same as those of the cosmic microwave background (CMB) dipole anisotropy 
due to the growth of cosmic structures. In other words exploring the aberration effect may give 
us a clue to the horizon-scale physics perhaps related to the cosmic acceleration. In this paper we 
develop a method to explore the dipole angular modulation from the pixelized galaxy data on the sky 
properly taking into account the covariances due to the shot noise and the intrinsic galaxy clustering 
contamination as well as the partial sky coverage. We applied the method to the galaxy catalogs 
constructed from the Sloan Digital Sky Survey (SDSS) Data Release 6 data. After constructing 
the four galaxy catalogs that are different in the ranges of magnitudes and photometric redshifts to 
study possible systematics, we found that the most robust sample against systematics indicates no 
dipole anisotropy in the galaxy distribution. This finding is consistent with the expectation from 
the concordance A-dominated cold dark matter model. Finally we argue that an almost full-sky 
galaxy survey such as LSST may allow for a significant detection of the aberration effect of the 
CMB dipole having the precision of constraining the angular direction to ~ 20 degrees in radius. 
Assuming a hypothetical LSST galaxy survey, we find that this method can confirm or reject the 
result implied from a stacked analysis of the kinetic Sunyaev-Zel'dovich effect of X-ray luminous 
clusters in Kashlinsky et al. (2008,2009) if the implied cosmic bulk fiow is not extended out to the 
horizon. 

PACS numbers: 98.65.Dx 



The amplitude of the cosmic microwave background (CMB) dipole anisotropy is about two orders of magnitudes 
greater than the anisotropics at higher multipoles that are primarily generated during cosmic epochs until the last 
scattering surface of redshift z ~ 1100. It is widely believed that the dipole anisotropy is produced by the Doppler 
effect due to the relative motion between the Earth, i.e. an observer, and the frame where the CMB looks nearly 
isotropic (hereafter we call it the CMB rest- frame). The measured CMB dipole amplitude tells that the relative 
velocity has an amplitude of vcym— 370 km s^^ (1.23 x 10"'^ in the unit c = 1) [H, The peculiar velocity consists 
of the five vector contributions the motion of the Earth around the solar system barycenter (~ 30 km s^^), 
the motion of the solar system with respect to the Local Standard of Rest (LSR) the (hypothetical circular) 
motion of the LSR around the Milky Way (~220 km s"^ (the lAU 1985 recommended value), or ~250 km s'^Q), 
the motion of the Milky Way in the Local Group, and the motion of the Local Group with respect to the CMB 
rest-frame. The origin of the fifth component, the peculiar velocity of the Local Group, is still uncertain and has 
been under discussion over the past two decades (e.g., [1, and also see and references therein). This peculiar 

velocity is believed to be generated by the spatial inhomogeneities of mass (mainly dark matter) distribution in nearby 
large scale structures via gravitational instability as predicted in the cold dark matter (CDM) dominated structure 
formation scenario. Therefore the peculiar velocity field is expected to reflect properties of structure formation in the 
low-redshift Universe. For example, the peculiar velocity of the Local Group with respect to the CMB rest-frame is 
estimated to be ~ 600 km s~^, based on the result of 0]. This is greater than the rms amplitude of the peculiar 
velocity, ^ 470 km s~^, predicted from the linear theory of the concordance ACDM model. However, due to the 
difficulties in inferring the mass distribution from the observ ed g alaxy distribution, i.e. the galaxy bias uncertainty, 
the origin of the peculiar velocity is not yet fully understood [l0l - [l4| . 




I. INTRODUCTION 
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Then a naive question may arise; how special is the CMB rest-frame? One may think that cosmic large-scale 
structures have formed via gravitational instability over 13. 7G years in the CMB rest-frame. In other words, the 
Universe may still stay in the CMB rest-frame over the cosmic age [s^. However, a possibility has been discussed 
that the CMB rest-frame may not coincide with the rest frame of galaxy distribution due to super-horizon scale physics 
[isl [l6j. In fact, the present-day Universe is in the mysterious phase, the cosmic accelerating phase. The origin of 
the cosmic acceleration is one of the most profound problems in modern cosmology and physics, and it is sometimes 
discussed that the cosmic acceleration is perhaps related to the horizon-scale physics. For example, while dark energy 
is one possible explanation of the cosmic acceleration, dark energy should have spatial perturbations inevitably on 
horizon scales, if it is not a cosmological constant (e.g. [13, [13 ^^'^ references therein). Or the primordial power 
spectrum, produced in the inflationary era, may have weird behaviors on horizon scales such as the truncated power 
spectrum as speculated from the low CMB quadrupole amplitudes (e.g., (l9l- [21j '). Thus exploring the rest frame of 
the present-day Universe may give a clue to the horizon-scale physics (also see [l6[). Furthermore, this may give an 
independent test on the cosmological principle, the isotropy and homogeneity of the Universe. 

The rest frame of the present-day Universe may be defined by the frame where the galaxy distribution looks isotropic 
on a sufficiently large scale (hereafter we will call it the matter rest-frame or the rest- frame of the local Universe) . The 
relative motion of the Earth to the matter rest-frame causes a dipolc anisotropy modulation in the observed galaxy 
distribution [l^. The dipole anisotropy is caused by the two effects. First, photons from galaxies ahead of/behind 
the Earth are blue-/red-shifted by the Doppler effect, causing their fluxes to be brightened/dimmed and therefore the 
galaxies to be included/excluded in the magnitude limited sample. Secondly, special relativity predicts the aberration 
of angles causing the surface number density of galaxies to be enhanced/suppressed in the direction ahead of/behind 
us, even if the intrinsic galaxy distribution is perfectly homogeneous on the sky. If the motion has a similar amplitude 
inferred from the CMB dipole, the induced dipole modulation is small at a sub-percent level. Hence a large-area 
survey of distant galaxies or radio sources is suited for exploring this dipole anisotropy, because the intrinsic galaxy 
clustering that has greater amplitudes at low redshifts may cause a sig nificant contamination. There have been many 
attempts made to explore this dipole anisotropy from such surveys [22142^ . In particular Blake and Wall [IJ] analyzed 
the radio source distribution based on the NRAO VLA Sky Survey (NVSS) data |25|, and then claimed a possible 
detection of the dipole anisotropy that are consistent with the CMB dipole in the amplitude and direction within 2a 
and la levels, respectively. 

The purpose of this paper is to explore the cosmological dipole sig nature from the galaxy catalogs constructed from 
the Sloan Digital Sky Survey (SDSS) Data Release 6 (DR6) data [2g. The use of the photometric SDSS galaxy sample 
has several advantages: the sky coverage is large (about 20% of the full sky), the photometry is well-calibrated and the 
photometric redshift information of each galaxy is available. In doing this, we develop a method to explore the dipole 
modulation in a pixelized galaxy distribution properly taking into account the covariances due to the Poisson shot noise 
and the intrinsic clustering contamination as well as the partial sky coverage. In particular the photometric redshift 
information of SDSS galaxies is useful to reduce the clustering contamination from nearby structures. Furthermore, 
we will discuss how a planned large-area galaxy survey such as the Large Synoptic Survey Telescope (LSST) can be 
useful to explore the dipole anisotropy. 

Here it may be worth mentioning advantages and shortcomings of our method compared to the peculiar velocity 
field studies [^43 ■ This paper concerns with the dipole pattern in the galaxy distribution due to the observer's motion 
with respect to the matter rest-frame. Comparing this and CMB dipole anisotropics can address if the inferred matter 
rest-frame agrees with the CMB rest-frame within the measurement errors. This paper does not discuss the origin 
of our motion (or the motion of the Local Group) nor which structures in the Local Group cause our motion with 
respect to the CMB rest-frame, the so-called convergence depth. This question is one of the main questions discussed 
in the recent peculiar velocity field studies [ol-fl^. 

The structure of this paper is as follows. After reviewing the aberration effect on the galaxy counts of a given 
magnitude limit in Sec. |TT1 we will develop a method to explore the induced dipole modulation from a galaxy catalog 
properly taking into account the covariance and the partial sky coverage in Sec. IIIIl In Sec. IIVI we make an estimate 
on the detectability of the aberration effect from the galaxy distribution for hypothetical galaxy surveys, SDSS- and 
LSST-type surveys, assuming that the Earth's relative motion to the matter rest-frame has the same amplitude 
with the CMB dipole. Sec. |V] shows the main results of this paper. After defining the galaxy catalogs based on the 
magnitude range and the photometric redshift information and then estimating the covariance matrix of the pixelized 
galaxy counts on the sky, we will show the results for an exploration of the dipole anisotropy from the SDSS DR6 
galaxy catalogs where the amplitude and angular direction of the aberration effect are treated as free parameters. In 
Sec. |VT] we will also discuss the forecast for an LSST-type survey. Sec. IVIII is devoted to summary and discussion. 
Throughout this paper we will employ the concordance ACDM model that is specified by {h, J7a, ^m.h'^, ^bh'^, erg, n^) = 
(0.73, 0.762, 0.127, 0.0223, 0.74, 0.951) [13 . 
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II. EFFECT OF THE EARTH'S PECULIAR MOTION ON THE ANGULAR NUMBER DENSITY 

FIELDS OF GALAXIES 

Even if there exists the matter rest-frame where the intrinsic galaxy distribution looks perfectly isotropic, the 
peculiar motion of the Earth relative to the matter rest-frame induces an apparent angular modulation in the galaxy 



number density field on the sky [15|]. There are two effects that cause this modulation. The first is the aberration 
effect, causing the observed number density field to be increased or decreased in the angular direction forward or 
backward of the Earth's motion, respectively. The second is the Doppler effect: depending on the shape of spectral 
energy density of a galaxy, the Doppler effect causes the apparent magnitude of galaxy to be brighter or fainter, 
leading the galaxy to be included into or excluded from the magnitude limited sample. The net effect arises from the 
linear sum of these two effects in a case that the peculiar velocity is much smaller than the speed of light. 

In the following, we will in more detail model the effect of the Earth's peculiar motion on the angular number counts 
of galaxies defined for a given limiting magnitude. In so doing, we assume for clarity that the galaxy distribution 
is perfectly isotropic on the sky in the matter rest-frame, and in other words we ignore intrinsic anisotropics in 
the galaxy distribution arising from large-scale structure formation. We shall come back to the intrinsic anisotropy 
contamination later. 

Special Relativity predicts that, due to aberration of angle, the relative velocity of an observer to the matter 
rest-frame causes the angular position of a galaxy observed on the sky to be displaced: 



1 - sin sine/) 
tan(/)obs ——5— « 1— J' 1) 

cos (p + p COS (p + p 

where (j)ohs denotes the angle between the direction of the galaxy seen by the observer and the direction of the 
observer's velocity v, (p is the corresponding angle in the matter rest-frame, and /? = |v|/c. In the second equality on 
the r.h.s. we assumed /3 <C 1, and ignored the term of 0{/3^). In the following we will similarly ignore the contribution 
of O{0^) for simplicity. The conservation of the number of galaxies tells that the observer sees an angular modulation 
in the galaxy distribution as a function of the angular direction d on the celestial sphere: 

= « n(l -f^2/3cosa) , (2) 

«0obs 

where 6 — [con 9 cos (p, cos sin Lp, sin 9) for the spherical coordinates (wc have employed the continuous field limit, 
which is a good approximation for cases of interest) and n is the intrinsic number counts per unit steradian. The 
angle a is defined by cos a = • v, where v denotes the angular direction of the Earth's peculiar velocity v on the 
celestial sphere and is fully specified by two parameters. 

Next let us consider the Doppler effect on the galaxy counting. Besides redshift due to the cosmic expansion, the 
peculiar velocity of an observer relative to the comoving rest frame of galaxy distribution causes a photon emitted 
from a galaxy to be redshifted or blueshifted depending on the angular direction of the galaxy relative to the peculiar 
velocity direction. The observed frequency I'obs is related to the rest frame frequency v-^^at via 

Vohs l+l3cos(j) 

where the angle (f) is defined in the same way as in Eq. ([T|) . If we assume that an intrinsic flux density of a galaxy [in 
the units of erg s~^ cm~^ Hz~^] is simply given by a power law as S-ccstiy) oc v^, the conservation of photon number 
tells that the observed flux density Sohs{v) is expressed as 

•^obsi^'obsj = •^rcsti^'rcst j "1 "37 = V "Jrcst (I'obs j • 1,4 j 

^^rcst ^^^obs ^Cobs 

The apparent magnitude of the galaxy in a given fllter is obtained by integrating the flux density over a range 
of transmission frequencies of the fllter, and therefore the observed magnitude, mobs, is related to the rest frame 
magnitude, rrirest, via 

rriohs = m^cst - 2.5(1 - p) log^p -q. (5) 

If we assume that the intrinsic number counts of galaxies, where galaxies brighter than a given limiting magnitude 
miini are included in the sample, is simply given by 

n{m < mii,„) cx 10^™'"", (6) 
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with X being a numerical coefficient of order unity |28| . the observed number counts are found from Eqs. ([3]) and ([5]) 
to be 

n(0; m < mum) = n{m < miost.iim = mim + 2.5(1 - p) log^g ij) 
= r;2-5-(i-f)n(™ < m,,^) 

~ [1 + 2. 5a;(l — p)/3 cos a] fi(m < miini), (7) 

where the angle a is defined in the same way as in Eq. 

Hence, taking into account both the aberration effect ([2]) and the Dopplcr effect ^ up to the first order of /?, the 
observed angular number density field of galaxies for the limiting magnitude is expressed as 

n{9; m < mii,„) = n{m < rniim) 1 + 2/3 cos a , (8) 

with the modified /3 parameter defined as [G^I 

/3 = [l + 1.25x(l-p)]/3. (9) 

Note again that the dependence of n{6) on angular direction 9 comes through the relation cos a = 9 ■ ir. Eq. ([8]) 
shows that the Earth's motion relative to the matter rest-frame induces a characteristic dipole pattern in the angular 
number density field of galaxies. The dipole pattern is fully specified by 3 parameters: /3 and 2 parameters for v. The 
average of Eq. ^ over the whole sky indeed satisfies conservation of the total number of galaxies: 

— / dQa n{9) = n. (10) 

If the matter rest-frame is the same as the CMB rest-frame with respect to which the Earth is moving with the 
velocity P « (1.231 ± 0.008) x 10"^ in the direction (Z,6) = (264?31 ± 0?04 ± 0?16,48?05 ± 0°02 ± 0°09) in the 
Galactic coordinates [ll,[3, the dipole amplitude is 0{P) - 0(10-^) as the prefactor in Eq. ([9]) in front of /? is of the 
order of unity for the number counts of galaxies in optical passbands [28h30| . To detect the dipole pattern in the 
galaxy distribution, a galaxy survey with full sky coverage is ideally needed, while a partial-sky survey such as SDSS 
makes it less straightforward to explore the dipole pattern, because the galaxy distribution displays a smaller angular 
modulation over the sky region observed, as will be discussed below. 



III. METHODOLOGY: A TEST FOR DETECTING THE ABERRATION EFFECT 

In reality, a galaxy distribution actually seen is quite far from homogeneous, rather displaying rich, hierarchical 
structures on various distance scales - cosmic large-scale structures. To measure the aberration effect due to the Earth's 
peculiar velocity, we need to discriminate the effect from the inhomogcneitics due to large-scale structures. In addition 
wc have to take into account observational effects such as the survey geometry and the shot noise contamination due 
to a finite number of galaxies. In this section we develop a methodology for measuring the aberration effect from a 
wide-field galaxy survey, which will be applied to the SDSS data in subsequent sections. 

For an actual galaxy survey, the density perturbation field has to be estimated from the discrete distribution of 
galaxies. Since the dipole modulation of interest appears over an angular scale of tt radian and we are not interested 
in small angular scales, it is convenient to consider a pixelizcd map of galaxy distribution, where the pixel size greater 
than degree scales would be sufficient for our purpose. We define the number density of galaxies in the i-th pixel as 

nobs(0,)= V d^ngW^,-,{9,-9)Sl{9-9g,,), (11) 

where 9i denotes the angular position of the i-th pixel, given by 9i — (sin 9 i cos if i, sin 9 i simp i, cos 6 i) in the celestial 
polar coordinates, 6^(9) is the Delta function, 9g,j denotes an angular position of the j-th galaxy, and the summation 
runs over all the galaxies used in the analysis, j = 1, 2, . . . , A^gai.tot (.^gai.tot is the total number). Note that 
the integration J^^d^^lg is confined to the survey region with area D-s- The window function W(i'){9) defines the 
pixel shape normalized as J(P9 W(^i){9) = 1, where the subscript (i) is used for notational convenience to explicitly 
imply that the window function may change by pixel-to-pixel taking into account possible variations due to the pixel 
geometry, masking and so on. For example, a top-hat type window function is given by W(^i^{9) = l/(''''^pix(i)) if 
1^1 < ^pix(i)5 otherwise zero. 



5 



For our purpose wc need to deal with the density perturbation field defined as 

"-obs(^i) 



1, 



(12) 



Note that the iVgai,tot/f^s is an estimate on the average number density for a given survey. 

Accordingly we need to modify a modeling of the aberration effect such that the model prediction can be compared 
with the pixelized galaxy distribution above: the number density in the z-th pixel can be simply modeled from Eq. (|8]) 

as 



1 + 2/3 cos ai 



(13) 



where cos at = Oi ■ w. However, when a survey has a partial sky coverage, the conservation law ()10|) no longer holds. 
In other words, 



vis 



n{e) = n 1 + 2/3 



e 



cos V • I ^ n 



(14) 



where Vis is the survey area, and the integration range is confined to the survey region. Hence, the average number 
density n in Eq. p3l) needs to be estimated taking into account the partial sky coverage: the model density perturbation 
field, just like the derivation in Eq. ([T^ . is modified as 



(^modcK^'i) — 



1 + 2/3 COS ai 



(1/^^s) E;=T ^^pixO) (1 + 2/3 COS a, 

cos a,; - ^ X! ^Vi^]) 



- 1 



2/3 



(15) 



where the summation runs over all the pixels used (iVpix is the total number). Note that the denominator in the first 
line gives an estimate on the average number density for a case of the partial sky coverage, and the averaging of the 
aberration effect is weighted by the pixel area, ripix(j) • 

In this paper, to measure the aberration effect for a given galaxy survey, we employ a simplest statistic . From 
Eqs. ([T^ and p^ . the statistic can be given by 



x'(v,/3) 



E 



['5obs(0») 



□del 



(^.;v,/3) 



'5obs(0j 



□del 



(e,;v,/3) 



(16) 



where [C]^^ is the inverse of the covariance matrix as explained below. The best- fit model parameters for v and /?, 
3 parameters in total, can be estimated by minimizing the value with varying the model parameters freely, given 
the pixelized galaxy distribution. 

The statistical uncertainty in measuring the aberration effect is quantified by the covariance matrix C in Eq. (|16l) . 
Following [m (see Sec. 31; also see [l^), the covariance matrix is found to be given by 



(17) 



where 9ij = cos~^{9i ■ 9j), is the Kronecker delta function, A^g(i) is the number of galaxies contained in the z-th 
pixel, and Wg(9) is the angular two-point correlation function of the pixelized galaxy distribution. The first term gives 
the shot noise contamination due to a finite number of galaxies, while the second term quantifies the sampling variance 
originating from the intrinsic galaxy clustering in large-scale structure. Note that the second term is non-vanishing 
when i ^ j, which describes how the errors between different pixels are correlated with each other. The covariance 
matrix C has a dimension of A'pix x A^pix. 

The two-point correlation function Wg{6) in Eq. (jl7p can be expressed in terms of the angular power spectrum of 
galaxy distribution, Cg{l), as 



Wg{e,j, 



E 

1=1 
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(18) 
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SDSS-like (/,ky ^ 0.19) 


7838 


487 


~ 0.3 


5.3 X IQ-* 


1.3 X 10"^ 


2.8 X 10"" 


0.40 


I0W-2: full-sky (/aky = 1) 


41200 


487 


~ 0.3 


2.4 X IQ-^ 


2.8 X 10"^ 


6.6 X 10-'' 


0.84 


LSST-like (/eky ~ 0.46) 


19000 


9.7 X 10'' 


1.2 


1.2 X IQ-^ 


2.1 X 10"" 


3.2 X 10"^ 


5.5 



TABLE I: A rough estimates of the expected signal to noise ratios for three survey configurations. The "SDSS-like" survey 
assumes the survey geometry, galaxy redshift distribution and galaxy number density of the SDSS DR6 for survey parameters. 
The "low-z full-sky" survey assumes the full sky coverage, but other parameters are the same as the SDSS parameters. The 
"LSST-like" survey assumes the configuration of the future LSST which will be explained in Sec. IVII The notations are 
as follows, /sky = f2s/(47r): the sky coverage of the survey. A'^gai.tot/^s: the average number density of galaxies per unit 
square degrees. Zm'- the mean redshift of the galaxy distribution. Aab: the dipole amplitude coefficient of the aberration effect 
assuming the CMB dipole amplitude, = 1.231 x 10"^, where the survey geometry is also taken into account. Alss: 
the dipole amplitude coefficient expected from angular galaxy clustering in large-scale structures. Ap: the dipole amplitude 
coefficient expected from the Poisson noise due to a finite number density of galaxies. S/N: the estimated signal-to-noise (S/N) 
ratio for measuring the aberration effect. 



where W(j-)(a::) is the Fourier transform of the i-th pixel's window function W(^i){x), and Pi{x) is the l-th order Legendre 
polynomial function. The symbol * denotes the complex conjugate. Assuming a linear bias model between the matter 
and galaxy distributions, which is a good approximation on angular scales of interest, the angular power spectrum 
Cg{l) can be expressed in terms of the underlying linear power spectrum of mass perturbations as 

2 f°° 

Cg{i) = -bi / p^{k)\m\'k^dk, (19) 

~ Jo 

with 

/■oo 

Ii{k) = / D{z)ng{z)ji{kr{z))dz, (20) 
Jo 

where bg is the linear bias parameter of galaxies, r is the comoving angular-diameter distance, D{z) is the linear growth 
rate normalized as D{z = 0) = 1 today, ng{z) is the redshift distribution of galaxies normalized as n{z)dz = 1, 
ji{x) is the Z-th order spherical Bessel function, and P^ik) is the linear mass power spectrum today. The CDM 
based linear perturbation theory provides secure predictions for P^{k) as a function of cosmological parameters [e.g., 
see . Once the CDM model and the galaxy bias are assumed and the galaxy redshift distribution is estimated 
or known, one can make a secure estimate of the sample variance contribution to the covariance. We will use the 
FFTLog code developed in [sj] to compute the z-integration in Eq. ([20]) (also see [1^). An alternative method is to 
use the angular correlation function measured from the survey itself, however, it is generally difficult to obtain an 
accurate measurement of Wg{9) on large angular scales we are interested in. 



IV. A ROUGH ESTIMATE ON THE SIGNAL-TO-NOISE FOR A GALAXY SURVEY 



Before going to the results, it would be useful to have a rough estimate on the signal-to- noise ratio for measuring 
the aberration effect from a given galaxy survey. As discussed in Sec. Illli the measurement is contaminated by the 
shot noise of discrete galaxy distribution and the galaxy clustering contribution in large-scale structures. Here we 
simply compare the expected dipole amplitude due to the aberration effect with those of the contaminating effects, 
taking into account the survey geometry (see [s^ for the details). A detailed derivation of the equations used in this 
section is also given in the appendix |Al 

As partially discussed in Sec. lIIIi an incomplete sky coverage dilutes the dipole modulation in the galaxy distribution. 
For the case of a partial survey coverage, the dipole amplitude of the aberration effect is estimated as 




The quantity 5*11 takes into account the dilution effect due to the survey geometry and is found to be given by 

I I' 

Su'^ J2 IWIP-^'I', (22) 

m— — / rn' — — l' 
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with the window function WJ]}"^ being defined in [33| as 

wj]}"^' EE Liig Yr{e)Yf*{e)M{e), (23) 



where YJ™ is the spherical harmonic function and M{0) is the mask function to define the survey region: M{0) = 1 
if the angular position 9 is inside the survey region, otherwise zero. 

The dipole component arising from the clustering distribution of galaxies can be characterized by the angular power 
spectrum, weighted with the survey window function. Following the method developed in [sB], the dipole amplitude 
can be estimated as 



Alss = \/ — . E^^'C^aO, (24) 

\ l'=l 

where the prefactor A/3/47r is from the definition of the first-order spherical harmonic function, Yi ~ y^S/Air cosO. 
The partial sky coverage causes the power spectra of all the I'-th orders, Cg{l'), to contribute to the dipole amplitude. 
The dipole amplitude arising from the shot noise due to the discrete galaxy distribution is similarly estimated as 

^ \/A7r^ -^TT fig ^ ^"^ 

Note that l/fig is the inverse of the average angular number density of galaxies per steradian. 

Hence, from Eas. (|2T|) . ((24| and (p5)) . the signal-to- noise ratio for measuring the aberration effect may be estimated 

as 

^ ^^"^ (26) 



The S/N simply assesses the dipole amplitude of the aberration effect relative to those of the contaminating effects. 

Table U gives the estimates on S/N for three types of galaxy surveys assuming the CMB dipole amplitude for the 
aberration effect, assuming /? = /? = 1.231 x 10~^. First we consider a 'SDSS-like survey' that mimics the SDSS 
galaxy distribution analyzed in this paper, including the survey geometry, the total number of galaxies and the 
redshift distribution. In this case S/N ~ 0.4, smaller than unity, meaning that the aberration effect is difficult to 
measure. Comparing the dipole amplitudes, Aab, ^lsSj and Ap clarifies that Alss gives a significant contamination 
to an extraction of the aberration effect. To study the impact of the partial sky coverage, the row labeled as 'full-sky' 
shows the results when considering the full-sky coverage but keeping other survey parameters to be the same as those 
for the SDSS-like survey. The aberration signal becomes more significant, however, the clustering dipole amplitude 
contamination is still significant, resulting only in a slight improvement as S/N ~ 0.8. On the other hand, the LSST- 
like survey that probes the galaxy distribution up to much higher rcdshifts and with a wide sky coverage /sky ~ 0.46 
is found to allow for a significant detection of the aberration effect, S/N ~ 5.5, more than 5cr. Therefore these results 
show that it is more important to probe the galaxy distribution at higher rcdshifts for measuring the aberration 
effect mainly because of the following two reasons. First, the galaxy distribution is less evolving and more in the 
linear regime at higher rcdshifts, therefore the clustering dipole less contaminates to the dipole signal. Second, while 
the cosmic structures at smaller distance scales are more rapidly evolving towards the nonlinear regime in the CDM 
structure formation scenario, the small-scale structures at higher rcdshifts are viewed by an observer with smaller 
angles, less contributing to the dipole amplitude appearing at large angular scales. 



V. APPLICATION TO SDSS DR6 PHOTOMETRIC GALAXIES SAMPLE 

The SDSS 6th Data Release (DR6) [1^ covers about 8000 dcg^ of sky area and contains over 200 million of 
objects with photometry in five pass bands: u, g, r, i, and z (39l.l40|. The SDSS galaxy sample would be a most suitable 
data set to explore the aberration effect because of the well-calibrated, homogeneous photometric and astromctric 
properties. 



A. Constructing Galaxy Samples from the SDSS Photometric Galaxy Catalog 

We construct the galaxy sample from the SDSS photometric catalog following [4l|. There arc several uncertainties 
in the photometric calibration that may cause artificial angular modulations in the galaxy distribution, such as an 
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imperfect correetion of Galactic dust extinction and the contamination of stars to the sample caused by an imperfect 
start-galaxy separation. We will also use photometric redshift (hereafter, simply photo-z) information of each galaxy 
to define our galaxy sample. We will below briefly describe the definition of our galaxy sample. 



1. Dust extinction correction 



An inaccurate correction of Galactic dust extinction may cause artificial angular modulations in the galaxy distri- 
bution that is defined for a magnitude-limited sample. In the SDSS database, each object has information on not only 
its photometric properties but also the dust extinction estimated based on the Galactic dust extinction map in |42| 
(hereafter SFD), Ax,sfd, {x = u,g,r,i, and z). As carefully investigated in [4l| . there may remain a small systematic 
bias in the extinction correction, especially in the field with Ar^sFD < 0.1. However, the possible systematics has not 
been yet resolved, so we adopt the magnitudes where dust extinction was corrected for based on the SFD map. 



2. Star-galaxy separation 

A secure galaxy-star separation is also important, because star contamination to the galaxy sample likely causes 
angular modulations in the galaxy number counts towards the Galactic plane. We constructed a photometric galaxy 
sample taking the following three steps: 

1. False objects were discarded using photometric processing flags. Namely, we removed the photometric objects 
that have saturated fluxes, were observed during bad sky conditions or arc identified as fast-moving objects. 

2. Masked regions were excluded, where the masked regions are defined from regions labeled as "BLEEDING" , 
"BRIGHT STAR" , "TRAIL" , or "HOLE" . 

3. The magnitude range to define the galaxy sample was employed to ensure a reliable star-galaxy separation. 



The details on Step 1 and 2 can be found in [41|. In Step 3 we employed the two i-band magnitude ranges: 19.1 < 
nii < 19.6 and 19.6 < nii < 20.1, respectively. 



3. Reducing a clustering dipole 

Nearby non-linear structures are viewed by an observer on relatively larger angular scales after projection. This 
causes a contamination to our seeking of the aberration effect because the galaxy clustering on small scales contami- 
nates to the dipole amplitude for a survey with partial sky coverage. For example, the apparent galaxy distribution 
in the nearby large scale structures, which extends up to ~ 60 — 200 Mpch~^ in radial distance, is away from the 
CMB dipole direction only by ^ 10 — 20 degrees on the sky [e.g. [l0l - [T2| . In addition, from the SDSS galaxy catalog 
itself, the apparent galaxy concentration at z ^ 0.08, the so-called "Sloan Great Wall", has been found [4j], and is 
away from the CMB dipole direction only by ^ 10 degrees. Thus including such nearby non-linear structures in the 
galaxy sample may apparently enhance large-angle clustering amplitudes in the galaxy distribution toward the local 
structures, which in turn prevents us from detecting the aberration effect due to the Earth's peculiar motion relative 
to the rest-frame of the SDSS galaxy distribution at typical redshifts z ~ 0.3. Hence, to maximize a chance to detect 
the aberration effect on the SDSS galaxy sample, it is desirable to remove such nonlinear structures at low redshifts 
as much as possible. 

For the reasons mentioned above, we use photometric redshift information to define a secure galaxy sample. Fig. [T] 
shows the photometric redshift distribution of SDSS galaxies per unit stcradian as a function of i-hand magnitudes. 
One can see that the SDSS photometric galaxy sample has a typical redshift of z ~ 0.3. However, there also appears 
too large population of galaxies at z ^ 0.05: since the redshift range covers a very small volume, the galaxy population 
is very likely to be largely contaminated by the outliers of photometric redshift estimates. To avoid this contaminating 
population, we employ the lower cutoff Zph,min = 0.1 for our galaxy sample. Fig. [1] also shows that there are not many 
galaxies beyond z = 0.9. Furthermore, in order to study the impacts of low-redshift nonlinear structures and galaxy 
magnitude cut, we study the four galaxy catalogs listed in Table HIl Note that the magnitude cut difference in the 
bright and faint samples are intended to study a possible contamination of imperfect galaxy-star separation. 
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Galaxy Sample 



Bright-Shallow (BS) sample 
Bright-Deep (BD) sample 
Faint-Shallow (FS) sample 
Faint-Deep (FD) sample 



Magnitude range 



2ph range Ug [deg 
(NGH) 



Ug Ug 

(SGH) (NGH+SGH) 


ng(SGH)/n9(NGH) 


329.1 337.8 


0.972 


501.7 502.3 


0.999 


431.1 451.2 


0.951 


774.8 802.2 


0.963 



19.1 < mi < 19.6 
19.1 < mi < 19.6 
19.6 < mi < 20.1 
19.6 < mi < 20.1 



0.1 < 2ph < 0.4 
0.1 < Zph < 0.9 
0.1 < Zph < 0.4 
0.1 < Zph < 0.9 



338.6 
502.4 
453.1 
804.9 



TABLE II: The definition of the four galaxy samples used in the analysis of this paper. The columns show the following: 
the i-band magnitude range imposed to define the galaxy sample; the photometric redshift range imposed to define the galaxy 
sample; the average number density of galaxies for the Northern Galactic hemisphere (NGH) region with area 6928 deg^; the 
average number density for the Southern Galactic hemisphere (SGH) region with area 679 deg^; the average number density 
for the whole NGH-I-SGH region with area 7607 deg^; the ratio of the number densities of NGH and SGH regions. 



M 50 
T3 




19.1<m.;<19.6 
19.6<m,<20.1 



0.2 0.4 0.6 0.8 

photometric redshift: 



FIG. 1: The photometric redshift distribution (0 < 2: < 1) of galaxies used in this analysis, from the SDSS DR6 data. The 
solid curve shows the distribution for galaxies with i-band magnitude in the range 19.1 < rrii < 19.6, while the dotted curve 
for 19.6 <mi< 20.1. 



4- Constructing the galaxy number count field 

Our method described in Sec. IIIII is applicable to the pixelized data of galaxy distribution. For convenience, as 
done in [4l| . we employ the same pixelization as that used in the SFD dust extinction map, which is given in the 
format of pairs of 4096x4096 pixel Lambert projections, for each of the Northern and Southern Galactic hemispheres 
(hereafter NGH and SGH, respectively). Note that the pixel size is (2'.372)2 (sec Appendix C in SFD). The SDSS 
galaxy distribution has a partial sky coverage, so wc do not use the pixels that arc not included in the survey region. 
Also for safety to avoid the pixelization effect, wc do not use the pixels that reside at boundaries of the SDSS survey 
regions. For the remaining pixels, dust extinction effect on galaxy magnitudes is corrected for based on the SFD 
map, and then the galaxy number counts is computed in each pixel for each galaxy sample. However, since we 
are interested in the dipole (large-angular) anisotropy of galaxy distribution, we do not need small-scale clustering 
information. Therefore for computational convenience we create a coarser pixelized map of the galaxy counts in which 
each pixel is defined by combining 73 x 73 square-shape neighboring pixels. The area of each pixel (without masking) 
is 8.34 deg^. Furthermore, the SDSS regions contain masked regions as described before. Properly taking into account 
the masking effect, each pixel is assigned to the area of unmasked region and the combined galaxy number counts. 
Here for safety we do not use the pixels whose effective area is smaller than 10% of the pixel area after masking. As 
a result, the galaxy number counts we will work on are given in 983 and 168 pixels for the NGH and SGH regions, 
which have areas of 6928 and 679 deg^, respectively (1148 pixels with 7607 deg^ in total). 

Table HIl summarizes the average number density for the four galaxy catalogs. As mentioned in Sec. IV A^ above, 
we constructed the bright and faint samples to see effects of incomplete star-galaxy separation, expecting the faint 
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FIG. 2: The pixelized galaxy distribution of the BD sample in the Northern Galactic hemisphere (NGH: left panel) and the 
Southern Galactic hemisphere (SGH: right), respectively, based on the Lambert projection. The gray-color scales correspond 
to the number densities of galaxies in each pixel per unit square degrees as indicated by the right-side bar. In each panel the 
outer light-gray circle represents the zero galactic latitude {b — 0°), while the inner circle b = 45°. The numbers labeled along 
the 6 = 0° circle show the galactic longitudes. 

samples suffer from more contamination. On the other hand, the shallow and deep samples are prepared to see effects 
of the large scale structure contamination, expecting the shallow samples have more contaminant. Other than these 
systematics, the average density of galaxies is found to be larger in the NGH region than in the SGH region. Although 
we are not sure for the significance due to the limited sky coverage of SDSS survey, the constraints on the aberration 
effect is found to be sensitive to an inclusion of the SGH region into the analysis as will be shown below, because the 
number density difference, if it is not real, mimics the dipole modulation of galaxy distribution. 

Fig. [2] shows the pixelized galaxy distribution of the BD sample in the NGH (left panel) and SGH (right), respec- 
tively. The survey geometry can be clearly seen: the NGH region has a much greater coverage than the SGH, and 
the SGH regions have three survey stripes. 



The x^-estimation of the aberration effect, given by Eq. (|16|) . requires an estimate of the covariance matrix, Cy , 
where the indices i and j run over the pixels of galaxy density map. 

As can be found from Eq. (jl7p the diagonal components, Ca, consist of two contributions: the Poisson noise arising 
due to discreteness of galaxy distribution and the sampling variance arising from the density fluctuations in large-scale 
structures. The Poisson noise contribution can be directly computed from the galaxy number counts in each pixel. 
On the other hand, the sample variance depends on the variances of galaxy density fluctuations of pixel scales, and 
the computation needs a few cares. Firstly, the areas of pixels are not uniform due to masking, and may vary from 
0.834 to 8.34 deg^ as described above. Secondly, a theoretical estimate of the sample variances (see Eqs. [T8land fT9)) 
involves several uncertainties: galaxy bias uncertainty and nonlinear clustering uncertainties corresponding to the 
pixel scales. In other words the linear mass fluctuations and the linear galaxy bias very likely break down on the 
relevant scales. Thus we instead estimate the sample variances directly from the SDSS galaxy catalog itself. To do 
this we used the original finest pixelized data of (2.372)^ arcmin^ in order to estimate the variances as a function of 
the smoothing scales, where the smoothing is done by combining neighboring pixels of square-shape region. Note that 
the smoothing angular scale is simply estimated from the area of combined pixels 

Fig. |3] shows the estimated variances for the four galaxy samples over the entire SDSS survey region, where the 
variances are estimated for 18 different smoothing scales. We then assign the sample variance of an arbitrary pixel 
area by spline-interpolating the 18 data points. Thus we ignore the hole effect or pixel geometry due to masking 
for simplicity. Since the galaxy sampling region within each pixel in the galaxy catalog is not necessarily connected 
due to masking, the sample variance estimated in our method arises from the density fluctuations of the smallest 



B. Computation of covariance matrix 
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FIG. 3: The variances of galaxy number density fluctuations as a function of the smoothing scales, for the four galaxy catalogs 
defined in Table HH The variances are estimated from the pixelized galaxy catalogs with varying pixel sizes (here 18 different-size 
pixels are computed), and the smoothing scale in the horizontal axis is simply estimated from the pixel area as ^sm — \/ ^^pix / 
(see text for the details). Note that, for illustration purpose, the results for the FS and BD samples are slightly shifted in the 
horizontal direction. The error bars around each data points are estimated based on the bootstrap resampling method. 



smoothing scales for a fixed area and would be greater than the aetual sample variance. Thus our estimate of the 
sample variance is somewhat conservative. 

Recall that, as given by Eq. (jl7p . the Poisson noise contribution to the covariance is given by l/(ngr2pix(i)) for a 
pixel with area ripix(i) {n is the average number density), where ripix(i) varies from 0.834 to 8.34 deg^. Table HI] tells 
that, for our galaxy catalogs, the Poisson noise is in the range of 0(10"'*) — O{10^^). On the other hand, since the 
pixel areas in the range of 0.834 to 8.34 deg^ correspond to the smoothing scales from 0.52 to 1.63 deg in Fig. [3l the 
sample variance is found to be O(10~^). Therefore, for our galaxy catalogs, the sample variance dominates over the 
Poisson noise in the diagonal covariance components, more then by a factor 10, i.e. (samp. vari.)>10 (Poisson noise), 
thanks to the enormous size of SDSS galaxy catalog. 

An estimate of the off-diagonal covariance components, Cij{i 7^ j), is more straightforward, where it contains 
only the sample covariance contribution. Since different pixels of our galaxy catalogs are separated by more than 
2.89(« Vs. 34) degrees, the cross-correlation between the galaxy density fluctuations of different pixels are safely 
considered to be in the linear regime. Therefore we use Eq. ([TSl) to estimate the off-diagonal covariances. The 
computation requires several ingredients. For redshift projection, we used the photometric redshift distribution of 
galaxies in Fig. [T] As for the linear mass power spectrum we employed the transfer function given in [4J] assuming 
the WMAP-3year cosmology [13]. As can be found from Eq. (|18p the pixel window function is needed to specify; 
we adopted a top-hat type window function W{6) = l/(7r0j[j) where 6th is simply estimated from the pixel area as 
9th = \/f^pix(i) [6l|. Finally we also need to specify the linear bias parameter of each galaxy sample, bg. We 
estimated bg by fitting the variances at large angular scales 1.5° < < 2.6° to the linear theory predictions, resulting 
in the linear bias parameters bg ~ 1.0 for both the BS and FS sample galaxies, while bg ~ 1.2 for the BD and FD 
samples. 



C. Expectation 



Before going to the results, the upper panel of Fig. 2] shows the expected dipole modulation pattern of galaxy 
distribution in the SDSS survey region, if the aberration effect on the SDSS galaxy distribution is in the same 
amplitudes and angular direction as in the CMB dipole, i.e. if the Earth's peculiar velocity to the SDSS galaxy 
distribution is /?cmb = 1-231 x 10""^ in the direction (^cmb,^cmb) = (264?31, 48°05) in the Galactic coordinates. 
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FIG. 4: Upper panel: An expected dipole pattern of the galaxies number density field for the SDSS survey region. Here we 
assumed the dipole amplitude and direction as those of the CMB dipole: the velocity amplitude /? — 1.231 x 10~^, and the 
angular direction is denoted by the black square symbol. We set /? = /3 in Eq. (|9| for simplicity. The color scales denote 
the number density fluctuations of galaxies as indicated by the right-hand color bar. Lower panel: The expected accuracy of 
constraining the dipole angular direction by the fitting (Eq. [16] or I27p , where the covariance matrix for the FD sample is 
used as a representative example. The contours denote the (Ax^) distribution over the sky, and show that the angular direction 
is not constrained: all the directions are within the la region ((Ax^) < 2.3 for two parameter case). The cross and the plus 
symbols indicate the minimum and the CMB dipole direction, respectively. The (Ax^) value corresponding to each contour 
line is indicated by the numbers on the line. 



For simplicity we ignore the Doppler effect, i.e., set /3 = /3 in Eq. ([9]). The plot shows the O(10^'')-level density 
modulation can be expected, however, the partial sky coverage of the SDSS region seems to significantly obscure the 
characteristic pattern and prevent the detection. It may also be worth noting that the CMB dipole direction predicts 
a negative number density contrast in the SGH region, which may be indicated from our galaxy catalogs as shown in 
Table nil although it may be just a coincidence. 



Assuming the dipole density modulation shown in the upper panel of Fig. 31 we then compute the mean values 
(minus the degrees of freedom, or "d.o.f.") with varying free parameters (v,/3). In other words, using Eq. p6|) we 
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FIG. 5: An expected marginalized likelihood of the dipole amplitude /3 is shown for the SDSS galaxy distribution depicted in 
Fig. [l] The SDSS catalog allows for an accuracy of /? = O(10~^) for the aberration effect search. 



compute 



^{Sassnmcd{Oi)) - 6niodc\iOi;'V, (3) [C]^j^ ((^assumed (^j )) - <5modcl (^j ; V 



where (/) denotes an ensemble average of the quantity / and 

(<5assumod(^j)) = 2/3cMBVcMB ' dj- 



(27) 



(28) 



(By subtracting the degrees of freedom, this quantity becomes zero when our free parameters (v, /3) coincide with the 
hypothetical observation parameters (vcmb,/3cmb) hi this case). The survey geometry of the SDSS DR6 are taken 
into account when computing 6j, the covariance matrix and i5modei- 

The lower panel of Fig. |4] then gives a more quantitative estimate, showing the accuracy of determining the 
dipole angular direction from the expected number density modulation. To be more explicit, this plot shows the 
theoretically expected mean difference between the best-fit model and models with varying model parameters to 
specify the angular direction of the peculiar velocity, marginalized over the velocity amplitude computed from Eq. 
(|27l) : (Ax^(Z, b)). We used the covariance matrix for the FD galaxy sample (see Sec. IV B|) in order to compute the 
Ax^. The figure clearly shows that (Ax^) is smaller than 1 over the entire sky. That is, the aberration effect is very 
difficult to detect, if the peculiar velocity to the SDSS galaxy is similar to the CMB dipole amplitude. In other words, 
the covariance. especially the intrinsic galaxy clustering contamination, is so significant compared to the aberration 
effect. This result is consistent with a rough estimate shown in Table HI Fig. [5] shows the marginalized probability 
distribution of the dipole amplitude parameter j3. Only an upper limit on /3 at the level /? ^ 10~^ is likely to be 
obtained from the SDSS catalog. 



D. Results 



Now let us move on to the measurement results. Here in the analysis on the real data, we first count number of 
galaxies in each pixel and compute the density perturbation field using Eq. (jl2p . Second, choosing model parameters 
(v, /3), we compute (Smodci using Eq. (fTSj) and construct x^(v, /3) using Eq. (fT6|) for those particular model parameters. 
Varying the model parameters, we iterate the second step, searching for the parameters that minimize the value. 
The value of f3 is varied in the range of < /3 < 7.0 x 10~^ and the values of (/, b) are varied over the entire sky 
directions. We used the covariance matrix computed based on the method in Sec. IVBI for each galaxy sample. 

Table IIIII summarizes the x^-fitting results for the four galaxy samples listed in Table |lll The second- and third- 
columns in the table give the minimum x^ value, Xmim for the best-fit model and its model parameters. The fact 
X^ ~ (pixel number) implies that our covariance matrix estimate is reasonable. The table shows that the best-fit 
values of /3 are of the order of O(10~^) as expected from Fig. [51 Table Hill also shows the x^-differences, Ax^, obtained 
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Galaxy Sample 



best-fit model: Errors: (A/3, A6I) Xcmb - Xn 



+0.0067 
-0.0054'°* ) 

+0.0059 inno\ 

-0.0057' ) 
/+0.0039 4no' 
I, -0.0034' 

(±0.0023, 33° 



2 2 
X)3 = ~ Xmin 



BS sample 

BD sample 

FS sample 
FD sample 



1080.3 (0.0127,310° ,30°) 

1211.4 (0.0087,290°, -10° 

1108.0 (0.0156,290° ,60°) 

1288.1 (0.0121,280°, 75°) 



5.85 

2.01 

18.7 
21.2 



7.22 

2.35 

21.7 
25.6 



TABLE III: The results for the fitting for the four SDSS galaxy samples defined in Table [III The second column labeled 
as "Xmin" denotes the minimum value for the best-fit model, and the third column gives the best-fit parameters. The 
fourth column shows, for each sample, the la error for /? and the la uncertainty of determining the dipole direction in radius 
from the best-fitting angular position. The fifth column gives the x^-difference between the best-fit model and the case that 
the peculiar motion induced dipole has the same amplitude and angular direction as those of the CMB dipole, given as 
(/3cMB, ^CMB, 6cmb) = (1.231 X 10~^, 264?31, 48?05). The sixth column shows the x^-difference compared to no aberration 
effect case, 13 — 0. 




FIG. 6: The likelihood probability distribution for the dipole amplitude parameter /3 of the aberration effect, marginalized 
over the angular direction parameters {l,b), for the four galaxy samples. The FS and FD samples which preferentially include 
fainter galaxies prefer a non-zero 13 with the amplitudes O(10~^). The other two samples are consistent with zero /3. 



by comparing the Xmin with two models: one assumes the CMB dipole amplitude and direction for the aberration 
effect, and the other is no aberration effect {(3 = 0). For both cases, the differences are modestly large, except for 
the BD sample, implying that the dipole angular modulation pattern, which is not necessarily the aberration effect, 
is marginally detected. 



Fig. |6]is the marginalized probability distribution of the dipole amplitude parameter /3 for the four galaxy samples. 
The FS and FD samples seem to prefer non-zero /3, however, the likely amplitude is O(10~^), implying that the 
peculiar velocity of the Earth relative to the SDSS galaxy distribution is w ~ 3000km s~^. Such a large peculiar 
velocity is difficult to explain in the currently concordance ACDM model, and therefore we believe this is likely due 
to either the intrinsic galaxy clustering contamination or unresolved systematic effects, or both, as also demonstrated 
in Figs, m and [S] Wc will come back to this issue later in this section. 
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FIG. 7: Estimated marginalized probability of the dipole amplitude /3, where only the data in the Northern Galactic hemisphere 
region is used. 



Fig. [S] shows how the A^^ value, marginalized over the dipole amplitude parameter /3, varies as a function of the 
peculiar velocity directions on the sky. It would be interesting to find that the method developed here can constrain 
the velocity directions. Again, for the FS and FD cases, the dipole directions are constrained at more than Scr level 
(Ax^ = 11.8). Also interestingly, the preferred dipole direction is close to the CMB dipole direction. If a possible 
detection of the galaxy dipole is due to the intrinsic clustering contamination, the inferred direction is not necessarily 
similar to the CMB direction since we have removed the contamination due to the nearby large scale structures using 
the photometric redshift information. 

Now let us discuss possible causes of the indicated dipoles in our samples other than the ED sample. We construct 
the four galaxy catalogs to study possible remaining systematic effects due to the star-galaxy separation and the 
intrinsic clustering contamination, as are defined in Table HIl As shown in Table Hill and Fig. [51 the "faint" samples, 
which are supposed to be more contaminated by an imperfect star-galaxy separation than for the "bright" samples, 
tend to prefer the dipole directions at higher galactic latitudes. This is counter-intuitive, because the faint samples 
may suffer from more star contamination towards lower galactic latitudes, which in turn causes the dipole direction 
in the lower latitudes due to the enhanced overdensity. Hence the imperfect star-galaxy separation seems not a main 
source of the systematic effects. On the other hand, if comparing the "deep" and "shallow" samples, the shallower 
sample is supposed to be more affected by the intrinsic galaxy clustering contamination from more evolving nonlinear 
structures at lower redshifts, if the photometric redshifts are well reliable. In fact Fig. \E\ shows that the dipole 
amplitudes tend slightly larger for shallower samples. 

Another concern arises from the partial sky coverage. In particular, as discussed in Table [Til the SDSS galaxy 
samples tend to have a smaller average number density in the Southern Galactic hemisphere than in the Northern 
Galactic hemisphere, which may cause an (apparent) dipole anisotropy. Fig. [71 shows the probability distribution of 
/? if only the NGH region is used in the analysis. Now all four samples show broader distributions, and all but FS 
sample are consistent with zero /?. Fig. [H shows the same plots as Fig. [3 but here again we have excluded the 
data in the Southern Galactic hemisphere from our four samples in the analysis. The Icr contours in the four plots 
become smaller, indicating larger dipole amplitudes as consistent with the finding on Fig. [71 above. We also note that 
the minimum values of all four samples (green crosses in Fig. [3]) more or less tend to move toward the Northern 
Galactic pole than in Figs. [51 Hence the number density differences between NGH and SGH does move the direction 
of the detected dipoles in all the four samples toward the CMB dipole direction, yet the origin of the difference itself is 
unknown. As far as direction is concerned, the aberration induced dipole expected from the CMB dipole does predict 
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FIG. 8: The Ax^ distribution as a function of the dipole direction parameters {I, b) on the sky, marginahzed over the dipole 
amphtude parameter /3. The upper-left, upper-right, lower-left and lower-right panels are the results for the BS, BD, FS, and 
FD samples, respectively. The thick (black) contour corresponds to the la uncertainty region (Ax^ ~ 2.3). The 'x' symbol 
denotes the best-fit direction, while the symbol denotes the CMB dipole direction, (/cmb,&cmb) = (264?31, 48°05). 

such a difference between NGH and SGH as is indicated in the upper panel of Fig. 21 A data having a wider sky 
coverage, such as LSST, is desired to settle down this issue. 

The reason for all but BD samples, when excluding the SGH data, show the minimum values around the Northern 
Galactic pole may be attributed to imperfect Galactic dust extinction correction, although it is not perfectly clear how 
Galactic dust can produce the difference in the BS and BD samples (or, in other words, how effects of Galactic dust 
differ among different redshifts). In this paper we have assumed that the dust extinction correction by the SFD dust 
map is perfect, but it is well-known in the literature (see, e.g., (4l| ) that there may remain some residual systematics 
in the map. In summary, since the BD sample with and without the data in the SGH, which is the most robust 
sample against possible systematics among our four samples, indicates zero dipole, we conclude that the aberration 
effect is not detected from the SDSS DR6 in our analysis. 



VI. PROSPECT FOR A FUTURE WIDE-AREA SURVEY: LSST 

As we have so far shown, a survey with large sky coverage and sufficiently deep redshift coverage is required to 
detect the aberration effect. The planned galaxy survey, the Large Synoptic Survey Telescope (LSST) [i^ see also 
http://www.lsst. org/|, will provide us with a most promising opportunity to explore the aberration effect. Here we 
estimate the expected accuracy. 

To estimate the forecasts, we need to specify survey parameters. According to [46| . the survey area is assumed to 
be 18863 deg^, and the survey geometry is restricted to the region —75° < 6 < 15°, excluding the galactic disk region 
of —15° < 6 < 15° (also see the figure below). For the galaxy redshift distribution, we simply assume the analytic 
form given by Eq. (4) in [43] : 

ng{z) cx ^^e-^/^o (29) 

with zg = 0.4, which has the mean redshift z = 1.2. The galaxy distribution is normalized so as to have the mean 
number density of galaxies of 50 per unit square arcminutes [47j . The assumed LSST galaxy redshift distribution is 
shown in Fig. [10] along with the ones for the SDSS samples shown in Fig. [1] 
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FIG. 9: The same as Fig. [8] but here for these plots, we have used the SDSS DR6 data in the Northern Galactic hemisphere 
only. 




FIG. 10: The expected redshift distribution n„(z) of the LSST survey (dash-dotted line). The mean angular number density 
of 50 per unit square arcminutes is assumed [471 ]. Also shown for comparison are the galaxy redshift distribution of the SDSS 
bright sample 19.1 < rrii < 19.6 (solid line) and of the SDSS faint sample 19.6 < rrii < 20.1 (dashed line). 
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FIG. 11: The expected dipole pattern of the galaxies number density field expected for a hypothetical LSST-like survey, 
assuming the CDM dipole amplitude and direction as in Fig. |4l The diamond symbol indicates the CMB dipole direction. 
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dipole amplitude: f3 x 1 o"' 

FIG. 12: The expected marginalized probability of the dipole amplitude parameter /3 for the aberration effect on the LSST 
galaxy distribution in Fig. 1111 



Fig. [TT] shows the expected dipole anisotropy pattern in the LSST galaxy distribution, assuming the CMB dipole 
amplitude and directions. (For simplicity, we ignore the Doppler effect, i.e., set /3 = /3 in Eq. ^ to produce the figures 
in this section.) It is clear that the wide sky-coverage can nicely capture the dipole pattern. 



To quantify the detectability of the aberration effect with the LSST, we need to model the sample variance con- 
tribution to the covariance matrix. We first considered the pixclizcd map with 213 pixels; each pixel has an area of 
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FIG. 13: The expected mean Ax'^ (i-e., (Ax^) given by Eq. ((27])) distribution for tlie LSST survey. The contours are stepped 
by (Ax^) = 10, while the thick contour denotes the la region. The la accuracy of determining the angular dipole direction is 
~ 20 degrees in radius. 

88 square degrees or ^ 10 degrees scale. The pixel size is sufficiently large, so we use the linear theory to compute 
the sample variance contributions. However, to reduce the intrinsic clustering contamination from structures at low 
rcdshifts, we only include galaxies at rcdshifts greater than z — I. Also we assume the linear bias parameter bg = 1 
for simplicity. 

Figs. [T^ and [T3] show the forecasts, which show the marginalized probability of the /? parameter estimation and 
the expected mean Ax^ distribution over the sky, respectively. From both the plots, it is clear that the LSST may 
allow for a significant detection of the aberration effect: the accuracy of the /? parameter estimation is at a level of 
ct(/3) = 10~^, while the dipole directions can be determined with the precision of 20 degrees in radius. If we can have 
a full-sky survey with similar depth to the LSST, the ultimate precision is 10 degrees in radius. 

VII. SUMMARY AND DISCUSSION 

In this paper we have used the SDSS DR6 galaxy catalog to explore the dipole anisotropy pattern in the angular 
galaxy distribution caused by the aberration effect due to the Earth's peculiar motion to the matter rest-frame. The 
SDSS DR6 catalog, which covers about 8000 square degrees, allows us to construct the currently most reliable galaxy 
catalog for our purpose because the data has well-calibrated, homogeneous and secure photometric and astrometric 
properties. 

After developing a method to explore the aberration effect signal (see Sec. IIIip . we estimated the detectability of 
the dipole signal for SDSS- and LSST-type surveys assuming the CMB dipole amplitude and the concordance ACDM 
model properly taking into account the contamination effects (Sec. IIV|) . We found that the intrinsic galaxy clustering 
at low rcdshifts gives a significant contamination. Recalling that for a CDM model the galaxy distribution has greater 
inhomogeneities at smaller scales, the small-scale galaxy distribution at lower rcdshifts is viewed by an observer with 
larger angles, and the strong inhomogeneities significantly obscure the aberration effect if the Earth's peculiar motion 
to the matter rest-frame is similar to that for the CMB dipole amplitude, i.e. /3 = O(10~^). In other words, the 
dipole anisotropy due to the intrinsic clustering can be greater than the aberration effect for a low-redshift galaxy 
catalog. The Poisson contamination due to the discrete galaxy distribution is smaller, but not negligible even with 
such a huge number of SDSS galaxies. Therefore our results imply that the previous report on a possible detection of 
the aberration effect may be due to these contaminating effects (see also |48j). 

To remove the contaminating effect of the low-redshift galaxies as much as possible, we used the photometric redshift 
information to define a secure catalog of galaxies. We considered the four galaxy catalogs to study possible remaining 
systematic effects due to the star-galaxy separation and the intrinsic clustering contamination, as are defined in Table 
im As shown in Table UlIl and Fig. HI the "faint" samples, which are supposed to be more contaminated by an imperfect 
star-galaxy separation than for the "bright" samples, tend to prefer the dipole directions at higher galactic latitudes. 
This is counter-intuitive, because the faint samples may suffer from more star contamination towards lower galactic 
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latitudes, which in turn causes the dipole direction in the lower latitudes due to the enhanced overdensity. Hence the 
imperfect star-galaxy separation seems not a main source of the systematic effects. We also compared the "deep" 
and "shallow" samples by using the photometric redshift information. As described above, the shallower sample is 
supposed to be more affected by the intrinsic galaxy clustering contamination from more evolving nonlinear structures 
at lower redshifts. However, the dipole directions for both the deep and shallow samples are found to be similar, which 
indeed prefers a direction around the CMB dipole direction within the 1-cr region, except for the "faint-deep" sample 
(see Fig. 13). 

Among the four galaxy catalogs, the "faint-deep" and "faint-shallow" samples are found to give a possible detection 
of dipole amplitude given as /3 ~ 10~^ corresponding to the peculiar velocity v ^ 3000km s~^ (see Fig. |6]). These 
values are significantly larger than the CMB dipole amplitude by a factor 10. The peculiar velocity amplitude of 
lO'^km s~^ is similar to the virial velocities of a massive cluster, however, we know our Galaxy is not in such a cluster. 
On the other hand, the large-scale (relative) bulk flow of 1000km s~^ is at 2 — Scr deviations from the typical peculiar 
velocity of 470km s~^ predicted by the standard ACDM structure formation model. Furthermore, the previous studies 
on the peculiar velocity field for the Local Group have not found such a large peculiar velocity. Therefore we believe 
that a possible indication on the aberration effect on the SDSS galaxy distribution is not significant, and may be due 
to residual contaminating effects. 

We argued that the dipole direction may be caused, at least partly, by the apparent systematic difference of galaxy 
number densities between the SDSS survey regions of the Northern and Southern Galactic hemispheres (NGH and 
SGH, respectively), where the SGH region has a smaller average number density than the NGH (see Table HI)) . We 
found that the number density differences between NGH and SGH does move the direction of the detected dipoles in 
all the four samples toward the CMB dipole direction, however, the origin of the difference itself is unknown. The 
two regions are separated, and the SGH region has a much smaller sky coverage (679 deg^) compared with the NGH 
region (6928 deg^). A data having a wider sky coverage is desired to settle if this difference is real or due to some 
systematics in our analysis. 

Another concern is on the Galactic dust extinction correction. Indeed, when excluding the SGH data, the preferred 
dipole directions of all samples more or less move toward the Northern Galactic pole. Hence we suspect that Galactic 
dust may cause one part of our detection of the dipoles in the three less robust samples other than the BD sample 
(another possible part is the number density difference between NGH and SGH). Although we have assumed in this 
paper that the SFD dust map [42| gives us a perfect dust extinction correction, Yahata et al. [4l| pointed our 
remaining systematics in low dust extinction regions in the map. Therefore an indication of the dipole anisotropy we 
found is not yet conclusive. 

We quantitatively showed that a survey with almost full-sky coverage and sufficient depth is ideally 
needed to explore the aberration effect on the galaxy distribution. Future surveys, Pan-Starrs (49l . see also 
|http://pan-starrs. ifa.hawaii.edu/public/ or LSST [4a |. are such a survey and offer a chance to explore the signal. 
In this paper we showed that LSST may allow for a significant detection of the aberration effect. Even for the effect 
similar to the CMB dipole amplitude, the angular direction of the aberration effect can be determined wi th p recision 
of 20 degrees in radius. If the cosmic bulk flow of order ^ 1000km s~^ implied in Kashlinsky et al. (2008) [50| extends 
out to the horizon and therefore we are at rest in the matter rest-frame, the LSST survey will provide us a good 
opportunity for either accepting or rejecting such a bulk flow. Thus LSST may offer a unique chance to constrain 
the horizon-scale perturbations where some exotic physics related to cosmic acceleration such as dark energy and 
super-void may play a role. 

There are other methods/data-sets that allow to explore the horizon-scale peculiar velocity field relative to the 
Earth's motion. For instance, using a homogeneous catalog of quasars, such as that of SDSS, is advantageous to 
reduce the intrinsic clustering contamination due to the higher redshift coverage. However a clean identification of 
quasars from stellar population from imaging data is always problematic. Also the Poisson noise can be significant 
due to a much smaller number density of quasars. 

Another interesting method is using the kinetic Sunyaev-Zel'dovich (SZ) effect that is caused by the peculiar 
velocities of ionized medium that scatters off CMB photons causing the secondary temperature fluctuations by the 
Doppler shift. For example, the kinetic SZ effect can be extracted by measuring the CMB fluctuations. The advantage 
of this method is it allows to directly measure the line-of-sight component of peculiar velocity at the cluster redshift. 
In fact there are several attempts to measure the kinetic SZ effect by stacking the CMB fluctuations in the sky regions 
of clusters in order to explore the excess in the temperature fluctuations i50h52i |. In particular, Kashlinsky et al. 
(2008) [13 reported a possible significant detection of the large-scale bulk flow of X-ray luminous clusters, implying 
the bulk flow of ~ 1000km s~^, apparently similar to that indicated from our analysis on our faint samples within 
the uncertainties, in the angular direction not far from the CMB dipole direction (also see [H, lU)- Based on this 
surprising result, it was speculated that the large-scale bulk flow may originate from the non-standard horizon-scale 
perturbations, for example the tilt across the observable Universe due to the pre-inflationary inhomogeneities [l6| . 
This result is, very interesting though, still under debate, and a further analysis will be needed by using a homogeneous 
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massive sample of clusters. Moreover, a high-angular-resolution and high-sensitive CMB measurement has a potential 
to extract the kinetic SZ effect due to the ionized intergalactic medium (more exactly known as the Ostriker-Vishniac 
effect [HH). Thus it would be worth exploring large-scale bulk flow of the Universe by combining various methods in 
order to explore the possible new horizon-scale physics, which is very difficult to explore by other means. 
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Appendix A: Dipole Amplitudes 

This appendix briefly explains a method of estimating the amplitudes of the expected dipole moments in the galaxy 
density fluctuation field from three contributions: the Poisson shot noise, the large scale structure, and the motion 
of us (the Earth) with respect to the CMB rest frame. This method is used in Sec. llVl to obtain a rough estimate of 
a signal to noise ratio for detecting a dipole modulation. We first consider the case of a hypothetical all-sky survey, 
then a partial-sky survey. We follow the method of [H, 113, HI] with which readers' may consult for further details. 

1. All-Sky Survey 

This section considers a hypothetical all-sky survey and we will start with our definitions of the dipole moments. 
For a galaxy number density field n(0), we compute a spherical harmonic expansion of n{6) as 

ar = / dngn{e)Yrie), (ai) 

Jns 

where the integral is over the survey area i^s and * denotes a complex conjugate. We shall use the convention of [H^] 
for the spherical harmonics Yi™'{9). The dipole vector components in the three dimensional Cartesian coordinates is 
then defined as D = {—\/2TZe{a\), V2X m(a] ), a?) where TZe{f) and Xm{f) stand for a real and an imaginary part of 
some quantity /, respectively (See, e.g., [58|). Accordingly, the magnitude of the dipole is 

1 

m— — 1 

where the angle brackets denote an ensemble average. Given a™, we can estimate the magnitude of the dipole moments 
of various contributions. 
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a. Aberration Induced Dipole 

If the aberration is the only cause of the dipole in the galaxy number density field, the field has the form given by 
Eq. ([5]). Wc define the dipole amplitude due to the aberration effect as 



An 

yrigAab- (A3) 



where ng is the average angular number density of galaxies per steradian. By evaluating Eqs. (|A1[) and (|A2[) using 
Eq. ([8]), we can find Aab = 2/3 for an all-sky survey (and this is the reason for factoring out ng-\/47r/3 in the right 
hand side of the above equation). 

b. Large Scale Structure Dipole 
The coefficient a™ from the large scale structure (LSS) may be estimated as 

(l<LSsP)--^^&^ / P^;.{k)\h{k)\^k^dk = nlCS)- (A4) 
We then define the amplitude of the dipole due to large scale structure Alss as 



I /47r 

(li^LSsP) - y'3Q(l)n,EE^— n.ALss. (A5) 

c. Poisson Shot Noise Dipole 
For the Poisson shot noise, (a™p) = fig. We define the amplitude of the dipole due to Poisson noise as 



, /47r 

(|i?p|2) = ^3^=^_fi^Ap. (A6) 
2. Partial-Sky Survey 

The incomplete sky coverage may be taken into account by a tensor W^™™ [13, Ull ■ It is defined as 

wiir' = <f dngYr{e)YiV''*{e)M{^i), (a?) 

J Ait 

where an angular 'mask' M(J7) is defined such that M{Q,) = 1 over the area of the sky observed and Af(f2) = 
elsewhere. This tensor satisfies for any I 



^ E E E 1^-? = ^. (A8) 

in— — I I'—O m' — —l' 

With the tensor WIP''^'' , the coefficient of the spherical harmonic expansion for the incomplete sky coverage becomes 

(i^rn = E E i^r'i'iwfi'). m 

l'=Om' = -l' 

where a™ is the coefficients for the all-sky survey. We then define the amplitudes of the dipole fiuctuation in the 
galaxy number density field in terms of c™ instead of a™. Introducing coefficients 5*;/' as 

I I' 

Sw = E E 1^/""?' (AlO) 
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we define the dipole amplitudes from various contributions for a partial-sky survey as 

Alss 



\ 



(All) 



/47r V 47rng ' 

and 



= ^J^f > (A12) 



^ab = 2/3^^, (A13) 

which are the equations used to obtain the rough estimates of the signal to noise ratio of the aberration induced dipole 
in Sec. JVl 

Finally, we point out two things. Firstly, for a full-sky survey, W[P"^ = SjfiS^j^, and hence Sw = {21 -\- l)Slf,. As 
a result, we recover = 2/3 for such a survey as is expected. Secondly, we note that in a practical computation, 
we cannot complete the I' summation in Eq. (jAlip . We have made the summation up to I' = 500 for the SDSS DR6 
survey geometry and I' = 200 for the expected LSST survey geometry considered in this paper. Indeed, it is found 
that Eq. (|A8p satisfies to a good accuracy for these upper limit values for the ^'-summation. 

-. m— 1 — 500 m—l' ^ 

^ E E E l^i"™'!' = 0.9455^, 

m— — 1 I'—O m' ——I' 



Tn=l /'=200 Tn=i' „ 

= 0.9966- ^^^'^ 



Tn=l / =200 m=l 

^ E E E 1^1™ 



3 ^ ^ ^ ' " ' 47r 

m=— 1 ;'=0 m' =—l' 
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In General Relativity terms, the comoving hypersurface, where the energy flux of matter comoving with the cosmic 
expansion vanishes, has stayed the same from the last scattering surface to the present-day Universe. 
[60] In fact, we will consider flux limited as well as redshift limited samples in this paper (see Table HI)) . Because redshifts of 
galaxies ahead of our motion have systematically smaller than those behind us, there arises an additional dipolar variation 
in the number density. Up to 0(/?) and assuming we study galaxies in the redshift range zi < z < Z2, we can take account 
of this effect by writing 

/3 = (l-f 1.25x(l -p) + F{zi,Z2)) p. 

instead of Eq. Q where 

r 

(1 + 21)713(21) - (1 -I- 22)729(22) + / ng{z)dz 



2F(2l,22) 

where ng{z) is the galaxy redshift distribution in the matter rest frame. Assuming Eq. (|29|) for ng{z) with 20 =0.1 for 
SDSS or 20 = 0.4 for LSST, one can show that F(0. 1,0.4) = 0.33 (SDSS shallow sample), F(0. 1,0.9) = 1.4 (SDSS deep 
sample), or _F(l,co) — 0.91 (LSST). Hence F{z\,Z2) is positive and 0(1) in the case studies in this paper. 
[61] Here we simply ignore the square pixel shape and also the masking effect within one pixel. However, our approximation is 
good enough for the large-angle cross-correlations. 



